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We review the application of dynamical mean- field theory to Josephson junctions and 
study how to maximize the characteristic voltage I c R n which determines the width of a 
rapid single flux quantum pulse, and thereby the operating speed in digital electronics. 
We study a wide class of junctions ranging from SNS, SCmS (where Cm stands for 
correlated metal), SINIS (where the insulating layer is formed from a screened dipole 
layer), and SNSNS structures. Our review is focused on a survey of the physical results; 
the formalism has been developed elsewhere. 
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1. Introduction 

The DC Josephson effectfl, of a supercurrcnt flowing at zero voltage through a 
superconductor-barrier-superconductor sandwich, is one of the most fascinating 
macroscopic quantum-mechanicaLcffects in condensed matter physics. The orig- 
inal theoretical^ and experimental] work concentrated on superconductor-insulator- 
superconductor (SIS) tunnel junctions which havea hysteretic (double- valued) current- 
voltage characteristic. Ambegaokar and Baratoffl showed that for thin tunnel junc- 
tions, the characteristic voltage (product of the critical current at zero voltage I c 
and the normal state resistance R n ) is determined solely from the size of the super- 
conducting gap A; i.e. it is independent of the properties of the insulator. It was 
soon realized that the Josephson effect also occurred when the barrier was a normal 
metal (N) through the proximity effect as expressed by Andreev bound statesu. 
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Boguliubov and de Gennesu and Gor'kovO developed real-space formulations of the 
Bardeen-Cooper-ShriefferQ theory, that allow one to microscopically model inho- 
mogeneous Josephson junctions. The quasiclassical approach then progressed with 
the introduction of the (general) EilcnbergerO and (dirty) Usadet] equations, which 
were simplifications appropriate for ballistic and diffusive junctions, respectively. 
The challenge was to determine the proper boundary conditions for the partial 
differential equations that were appropriate for the experimental situations to be 
analyzed and to include self-consistency in the calculations. On the experimental 
side, new efforts focusing on ScS junctions^ (where c denotes a geometrical con- 
striction) were shown to be particularly suited to the quasiclassical approach, as the 
superconductivity was not diminished in the superconductor as one approached the 
interface of the geometrical constriction (and both the boundary conditions and t , 
self-consistency became trivial) . IBM embarked on a significant application effort 
where they created digital electronic circuits out of Pb and PbO based tunnel junc- 
tions and so-called latching technology (where the switching occurs between the 
zero voltage and finite- voltage parts of the I — V curve) . It became apparent that 
this latch-based technology would never be fast enough to be competitive with the 
limits of semiconductor technology, so the IBM project waned in the early 1980's. 

This was unfortunate, because a number of new breakthroughs occurred Jm 
the 1980's including the Nb-Al-AlO^-Nb process, developed at Bell Laboratoricsll3, 
which is used todayjjpr state-of-the-art low-T c based digital electronics; the Blonder- 
Tinkham-Klapwijkll3 model, which illustrated how multiple Andreev reflections can 
explain subharmonic gap structures in the I — V characteristic of SNS junctions; 
the development of rapid-single-flux-quantum (RSFQ) logico which showed how 
to make digital circuits run at the fastest possible speeds that superconducting 
electronics are capable of (and requires nonhysteretic single- valued I — V character- 
istics) ; and the discovery of high-temperature superconductorso which may be able 
to operate at speeds much in excess of a THz. Finally, the interest in mesoscopic 
superconductivitjJll and nanotechnologv have pushed efforts in the quasiclassical 
approach and Landauer-type approaches!!!! to the point where it is a well-developed 
tool to describe low-temperature junctions in the clean and dirty limits. Some 
milestones of the quasiclassical approach are the prediction that I C B«. can be in- 
creased over the Ambegaokar-Baratoff limit in point-contact junctionstSl, the ability 
to describe microscopic properties of SINIS junctions!*], the application pi multiple 
Andreev reflection theorycJ to submicron, self-shunted tunnel junctionsEil, and the 
understanding of how disordexjnodifies a clean junction through processes that take 
place at the Thouless energyE3. 

The success oflhe quasiclassical approach is dramatic and has been reviewed by 
numerous authorsEj. Nevertheless, there have been a number of new experimental 
results that examine Josephson junctions in regimes that lie outside the region of 
validity of the quasiclassical approach. These include (i) high-T c -based junctions^.! 
where the barrier is a correlated insulator with a charge redistribution occurring 
at the interface with the superconductor (grain boundaries, ion-damaged, interface 
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engineered, or Co-doped junctions); (ii) highly transparent SSmS junctions^ (Sm 
denotes a heavily doped semiconductor playing the role of a phase-coherent N) where 
the barrier can have its properties tuned by engineering the doping of the semicon- 
ductor (which is usually chosen to be either Si or InAs); (hi) SCmS junctionsc3 (Cm 
denotes a correlated metal) where the barrier is a correlated metal (or insulator) 
that lies close to the metal- insulator transition (such as TaN^); and (iv) mesoscopic 
short ballistic proximity-effect junctionstll where the transport through the barrier 
is ballistic, but the barrier thickness is less than the normal metal coherence length. 
These are junctions that may have technological importance and which provide a 
number of theoretical challenges. 

Our approach is complementary to the quasiclassical methods. We begin from a 
self-consistent many-body physics technique that automatically accounts for correla- 
tion effects and is best suited for short-coherence length superconductors (since thexe 
is no averaging over the Fermi length Xf)- It uses the dynamical mean held theorja 
approach, as modified for inhomogeneous systems by Potthoff and NoltingEa. We 
review briefly how these calculations are performed, but concentrate mainly on il- 
lustrating and ramlaining the resultsas applied to SNS and SNSNS structuresElj, 
SCmS junctionsEM3, SINIS junctional where the insulating layers are created by a 
mismatch of the Fermi energies of the S and N producing the ch arg e redistribution 
of a screened dipole layer, and ballistic proximity-effect junctionso. 

The coherence lengths in technologically interesting low-T^jjiaterials range from 
approximatelyEJ 40 nm in Nb, toEJ 24 nm in NbTiN films, tcEJ 5 run in Nb 3 Sn, to 
a rangecSof 1 — 5 nm for c-axis MgB2. In high-T c materials, the coherence length is 
typicallycj about 0.3 nm along the c-axis and 1—2 nm in the ab-plane. We choose 
to examine relatively short coherence length s-wave superconductors here, with the 
coherence length ranging from about 1 — 2 nm depending on the size of the unit 
cell (a more complete description is given below). 

We focus mainly on the figure-of-merit (or characteristic voltage) I c R n of the 
junction. Typical values in low-T c junctions lie in the range between 50 /zV and 
1.5 meV. The characteristic voltage determines the operating speed of a Josephson 
junction circuit element using RSFQ logic because the integral (over time) of the 
voltage pulse is equal to a flux quantum; hence the width of the pulse is inversely 
proportional to the height of the pulse which is determined by the characteristic 
voltage. This figure-of-merit has been analyzed in a number of different situations: 
(i) Ambegaokar and Baratofra showed that I c R n — 7rA/(2e) for tunnel junctions 
with vanishingly thin IJayers (A is the superconducting gap in the bulk); (ii) 
Kulik and OmelyanchukEJ showed that for clean point-contact SNS junctions one 
has I c R n = 7rA/e and for dirty (diffusive) point-contact junctions one has I c R n = 
27rA/(3e); Bardeen and Johnsorcd and IshiiH2l showed lhat long junctions have 
I c R n oc evp/L; and (iv) Freericks, Nikolic, and MillerEll showed that for clean 
wide three-dimensional SNS junctions, the charact eris tic voltage is limited by the 
product of the bulk critical current times the Sharvinc3 resistance, which typically is 
close to the above values. Experimentally, one finds a wide range of characteristic 
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voltages, especially as the barrier material is chosen to be a N, a Cm, a Sm, or 
a hybrid SINIS or SNSNS structure. One of the fundamental questions is how 
can one maximize the characteristic voltage while maintaining nonhysteretic I — V 
characteristics (tunnel junctions can be made nonhysteretic by shunting them with 
a low resistance metal shunt, which reduces the effective characteristic voltage, and 
thereby reduces performance, and requires more area for the shunted junctions on 
the chip). 

Another aspect of the problem that is important for hig h current density, short 
coherence length junctions is the issue of self consistency.OCj The superconducting 
order parameter A (or more correctly, the pair amplitude F which is the ratio of 
the order parameter to the interaction strength U; F — A/\U\) can be suppressed 
in the superconductor over a length scale on the order of the bulk superconducting 
coherence length (called the "inverse proximity effect"). In junctions where such a 
situation occurs, self consistency is critical for determining how the superconduc- 
tivity varies through the junction and in determining the size of the supercurrent 
when there is a phase gradient over the system. In addition, if the superconducting 
gap is large enough that the ratio of A to the Fermi energy /i (measured from the 
bottom of the bandlis nonnegligible, then there is intrinsic scattering induced by 
so-called A//x termsEJ. Our calculational methods are inherently self-consistent, as 
the DMFT requires self consistency to properly determine the effects of many-body 
correlations. This ensures current conservation within our junctions. 

2. Formalism 

The junctions that we model contain stacked infinite planes, consisting of super- 
conducting and barrier materials (B). The junction is divided into two pieces — (i) a 
bulk superconducting piece and (ii) a self-consistently determined superconductor- 
barrier-superconductor piece. This is illustrated schematically in Figure |l| where 
we show a self-consistent piece of 4 superconducting planes, 2 barrier planes, and 4 
superconducting planes, embedded in semi-infinite bulk superconductors on the left 
and the right. In our calculations, we always choose 30 self-consistently determined 
superconducting planes, and from 1 to 80 barrier planes (the word "barrier" is used 
to generically describe the "weak-link" material which can be a N, Sm, Cm, or a 
more complicated hybrid structure). The planar directions are denoted by x and 
y, while the inhomogeneous direction (perpendicular to the planes) is called the z 
direction. 

The system is discretized to a lattice (of lattice constant a) that represents the 
unit cells of the underlying ionic lattice of the junction^ This allows one to include 
arbitrary bandstructure and pairing symmetry effects.CJ We choose a simple cubic 
lattice with a nearest neighbor hopping integral t = 1 to describe the conduction 
band (we choose the hopping integral to be the same for the S and the B for simplic- 
ity; this choice can be relaxed at the cost of introducing additional parameters). The 
S is modeled by an attractive Hubbard model,c3 with an instantaneous attraction 
U = — 2. The Hubbard model is solved by an s-wave Hartree-Fock approximation, 
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Fig. 1. Inhomogeneous planar structure of the modeled Josephson junction. The junction is 
separated into two pieces: (i) a bulk superconductor piece and (ii) a self-consistently determined 
supcrconductor-barrier-superconductor sandwich. In calculations reported here, we always use 30 
S planes (about 8 times the bulk coherence length) and the barrier size ranges from 1 to 80 planes. 



which is equivalent to the BCS approximation^ except that the energy cutoff is 
determined by the electronic bandwidth rather than the phonon frequency.E-H The 
bulk superconducting transition temperature is T c = 0.11 and the superconducting 
gap at zero temperature is A = 0.198, which yields the expected BCS gap ratio 
2A/fcsT c w 3.6. The bulk coherence length £s is found by fitting the decay xif 
the pair amplitude F due to the "inverse proximity effect" at the SN interfaceHd 
We find £s ~ 3.7a, which agrees well with the BCS prediction of £s = ^v§,/(ttA) 
with vf, chosen as an appropriate average over the Fermi surface (for lattices with 
a lattice constant a « 0.3 — 0.6 nm, we find £s ~ 1 — 2 nm). The bulk critical 
current per unit area is / c ,bulk = 0.0289(2ei)/(?ia 2 ). The value of our bulk critical 
current density is slightly higher than the one determined by a Landau depairing 
velocity v c i = A/ftfep C? c ,buik = envd, where the density of electrons is n = k F /2Tr 2 , 
assuming a spherical Fermi surface) because of the possibility to have gaplesssu- 
perconductivity in three dimensions at superfluid velocities slightly exceeding E3 Vd 
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(note that Uf is direction-dependent for a cubic lattice at half- filling). Calculations 
on our junction are performed at a temperature of T = 0.01, which is effectively at 
the zero-temperature limit (T/T c w 0.09) for the superconducting properties. 

The barrier material will be chosen to be either a clean normal metal (which has 
U = 0) or an anneale d b inary disordered material described by the spin-one-half 
Falicov-Kimball modelJla The Falicov-Kimball model has two types of particles: (i) 
conduction electrons (which do not interact with themselves) and (ii) static ions, 
which can be thought of as classical particles that occupy a lattice site i if Wi = 1 and 
do not occupy a lattice site if Wi = 0. There is an interaction between the conduction 
electrons and the static ions denoted by Uf K , which results in an annealed binary 
distribution of the disorder. When the barrier is a clean normal metal, then one can 
calculate the coherence length of the junction to be £o = fivp /(ttA); in junctions 
where the hopping integral and Fermi level in the S and N are chosen to be the 
same, we always have £o = £s- This result does not hold in Sm barriers, because 
there is usually significant mismatch of the Fermi velocities, yielding different length 
scales. Nor does it hold in Cm junctions where the concept of a Fermi velocity is 
ill defined, and the length scale £o can only be estimated by mapping onto an 
effective disordered Fermi liquid and comparing the diffusive Thouless energy to 
the superconducting gap. 

The final interaction that we have in our system is a long-range classical Coulomb 
interaction, which generates a self-consistently determined electric potential Vf. 
This Coulomb interaction arises when the chemical potential of the S differs from 
that of the B. When these materials are assembled to form a Josephson junction, 
the overall chemical potential is fixed to be that of the S. The B layers then are 
at the wrong chemical potential, so a charge redistribution occurs, with charge 
shifting most near the SB interface. The value of the Coulomb potential V c at 
the crth plane is found by adding the potentials from every other plane that has a 
total electronic charge that differs from the bulk charge density of the respective 
plane (S or B). These potentials cause a local shift of the chemical potential equal 
to Vi = eVf , These potentials must be self-consistently calculated, so that the 
potential shift at a given plane, is equal to precisely the magnitude of shift needed 
to create the redistributed charge density of the plane (a more complete discussion 
will be given below). If the chemical potentials of the bulk S and B are equal to 
each other, then all Vi vanish. 

Hence the Hamiltonian of the junction is 

H = -J2 c t c jv + X] u * ( c 1t Ci T - ( c !i c a - 

+ E U i Kc Uio (wi - ~) + + AE F ) (cJ T c lT + c{ Cli ) , (1) 

ia ^ ^ i 

where ct (cj ff ) creates (destroys) an electron of spin cr at site i on a simple cubic 
lattice, Ui — —2 is the attractive Hubbard interaction for sites within the super- 
conducting planes, Uf K is the Falicov-Kimball interaction for planes within the 
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barrier, w; is a classical variable that equals 1 if a disorder ion occupies site i and is 
zero if no disorder ion occupies site i, Vi is the self-consistently determined Coulomb 
potential energy (if there is a charge redistribution), and AEp = Ep — Ep is the 
mismatch of Fermi levels in the S and B (AEp always vanishes in the S, but may be 
nonzero in the B). A chemical potential /i is employed to determine the filling. The 
superconductor is always chosen to be at half filling here (715 = 1), hence fi = 0. 

We employ a Nambu-Gor'kov formalismEffl to determine the many-body Green's 
functions. Details of the algorithm appear elsewhere. HEU'IE Since the junction is 
inhomogeneous in the z-direction only, we have translational symmetry in the planar 
directions. We begin by converting the three-dimensional problem into a quasi-one- 
dimensional problem, by using the method of Potthoff and Nolting.c3 We perform 
a Fourier transformation within each plane to determine the mixed-basis Green's 
function [defined in terms of two-dimensional momenta (k x and k y ) and the z- 
coordinatc of the plane] under the assumption that the electronic self energy is 
local, but can vary from plane to plane (in other words, the self energy has no k x of 
k y dependence, but does depend on z). For each momentum in the two-dimensional 
Brillouin zone, we have a one-dimensional problem with a sparse matrix, since the 
only coupling between planes is due to the hopping to each neighboring plane. The 
infinite "block-tridiagonal" matrix can be inverted by employing the renormalized 
perturbation expansion^, which calculates both the single plane and the nearest 
neighbor Green's functions. A final summation over the two-dimensional momenta 
produces the local Green's function and the Green's function for propagation from 
one plane to its neighboring plane. The DMFT is then employed to calculate the 
local self energy from the local Green's function and then the local Green's function 
is calculated from inverting the quasi-one-dimensional matrix. For the S or N, this 
amounts to just the Hartree-Fock approximation; for the Falicov-Kimball model, 
the exact solution corresponds to the coherent potential approximation (within the 
Nambu-Gor'kov formalism). These steps are repeated until the Green's functions 
have converged to a fixed point, where we have a self-consistent solution of the 
inhomogeneous problem that allows for nonuniform variations in both the pair-field 
correlations and in the phase. One important consistency check is total current 
conservation at each plane in the self-consistent region. There can be discontinuities 
in the current at the bulk-superconductor-self-consistent superconductor interface; 
the superconducting gap has always healed at this point, but there can be a jump in 
the phase (since this is far from the Josephson junction, it has a negligible effect on 
the results). This computational algorithm is a generalization of the conventional 
Boguliubov-de Gennes approach to allow for correlations within the barrier. 

This algorithm can be carried out for the normal state or for the superconducting 
state on the imaginary or real frequency axes. We work on the real axis in order 
to calculate the normal state resistance and the interacting DOS. Since we have a 
many-body system, we must use Knho's formula for the conductivity. Details for 
this calculation appeared elsewhereE9'El Our formalism calculates the conductivity 
by neglecting vertex corrections and evaluating the simple bubble diagram (which 
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is exact in the infinite-dimensional limit! 




3. Results for SNS junctions 

We begin by examining Josephson junctions where the barrier material is a normal 
metal, corresponding to U = 0, (w) = 1/2, and Ufk = 0.05. We choose the Falicov- 
Kimball interaction to be nonzero, because it makes the numerical computations 
converge more rapidly (by introducing damping), yet the transport through the 
normal metal remains ballistic, with a mean-free-path much longer than the barrier 
thickness. 

The first thing we will investigate is the critical current per unit cell. We can 
determine the dc Josephson current within an imaginary-axis calculation. We begin 
by calculating the supercurrent in the bulk superconductor when a uniform phase 
gradient is applied to the anomalous Green's functions (or equivalently the anoma- 
lous self energy). The bulk supercurrent increases monotonically with the phase 
gradient until the bulk critical current density is reached. Then we use this bulk 
solution for the boundary conditions of the superconductor and employ the compu- 
tational algorithm described above to calculate the local self energy (and Coulomb 
potentials, if relevant) on each plane. Next, we calculate the Green's functions (on 
the imaginary axis) that create an electron on the a + 1st plane and destroy an 
electron on the ctrth plane. The current passing through the ath plane is a simple 
integral of this Green's function! 3 !!. A strict convergence criterion is local current 
conservation — that the current through each of the self-consistent planes agrees to 
one part in 10 3 . Once the current is determined for a given bulk phase gradient, it is 
important to determine the total phase difference over the barrier, in order to calcu- 
late the current-phase relation. Since the barrier is spatially extended, and because 
the bulk superconductor has a uniform phase gradient over it, the total phase over 
the barrier must be carefully defined. We choose the following procedure: (i) first, 
we define the barrier region (of Nb planes) to correspond to the region that lies in 
between the midpoint of each SB interface (zl = N$ + 1/2, zr = N$ + Nb + 1/2); 
(ii) the phase difference over the barrier (of thickness L = Nsa) is determined by 
the total phase change over the barrier region, which corresponds to the uniform 
phase gradient LV0 plus a phase deviation term 54>{zp) — 64>(zl)', and (iii) since the 
phase is defined only at each lattice point, we define the phase deviation at zl and 
zr as the average of the phase deviation at the two neighboring lattice points { i.e., 
S(f>(zi J ) — {5(t)(Ns) + 5cf)(Ns + l)]/2 and similarly for the right hand side}. In most 
calculations, we can only determine the increasing part of the current-phase relation, 
where the current increases monotonically from to J c , but in some cases, we find 
that in the self-consistent region we can determine some of the "unstable branch" of 
the current-phase relation, where the current decreases as the total phase difference 
increases. In these cases, there must be a jump in the superconducting phase at 
the bulk S-self-consistent S boundaries (since the bulk supercurrent is always in- 
creasing with increasing bulk phase gradient). This does not affect our results much 
though, because the bulk S-self-consistent S boundary lies far from the barrier of 
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the Josephson junction, and the current is uniform within the self-consistent region. 

We plot the current-phase relation for a variety of SNS barrier thicknesses (1,2, 
4, 6, 8, 10, 14, 20, 30, 40, and 60) in Figure || (more correctly, we plot the current 
per unit cell area a 2 ). The dotted line is the bulk critical current density. One can 
see immediately that for thin junctions, where the I c approaches that of the bulk 
superconductor, the junction can never build too much total phase over it, so the 
current-phase relation is far from sinusoidal. As the total critical current of the 
junction decreases (as the thickness L increases), we find that the current- phase 
relation becomes more and more sinusoidal. This is illustrated further in the inset, 
where we plot the renormalized current-phase relation /(</>)/ I c - The deviation from 
sinusoidal behavior for high current junctions, with the critical phase difference 
being much less than tt/2, arises entirely from the self-consistencyll3. But it can 
be understood in a simple fashion, since the superconductivity is not too strongly 
reduced for a thin normal metal barrier, the majority of the phase difference is the 
bulk phase gradient multiplied by the thickness of the barrier. As the thickness 
gets small, the total phase across the junction must also be small, hence we get 
the behavior shown in Figure Note that the maximum of J (/b) always occurs 
below 7r/2. This difference from the analytical predictionaid'EJ of it for thick SNS 
junctions happens because of the self consistency.^ 

It is interesting to examine how the phase deviation behaves as the bulk phase 
gradient is increased from zero to the gradient corresponding to the critical current. 
Here we show such a plot for an intermediate thickness Nb = 20 in Figure |. These 
results are generic for most junctions. The phase deviation starts off negative, 
goes through zero at the center of the barrier, and then becomes positive (in a 
mirror image) of the negative function. What this tells us is that the majority 
of the phase change, or in other words the maximal phase gradient, occurs at the 
center of the barrier. This is where the superconductivity is the smallest. The 
critical current is determined by the maximal phase gradient that the central plane 
of the barrier can support and still maintain cur ren t conservation. Note further, 
the "notch" in the curves at the SN boundaries.^ This jump first smooths out, 
and then disappears as correlations are introduced. The phase-deviation curves are 
well-behaved and smooth for correlated SCmS junctions. In the case of thinner 
SNS junctions, the phase deviation curves develop an unusual "phase antidipole 
structure" for small current (where the slope of the phase deviation is negative at 
the center of the barrier), that disappears as the critical current is approachedll. 
The phase antidipole is. the opposite of the phase dipole behavior discovered in 
narrow SNS junctions. E2l We don't plot that behavior here. 

This behavior can be understood better if we directly examine the proximity 
effect for the Josephson junction. As described above, the pair-field amplitude F is 
a better measure of the superconductivity than the gap function A because the gap 
is always zero within the normal metal. The proximity effect is plotted in Figure ^ 
for a variety of junction thicknesses. The first thing to notice about the proximity 
effect is the suppression of F within the S as one approaches the SN boundary. The 
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Fig. 2. Current-phase relations for SNS junctions ranging from thin [Ng = 1) to thick (Ng = 60) 
barriers. The dotted line is the value of the critical current (for a unit square) in the bulk 
superconductor. Inset is the renormalized current-phase relation I((f>)/I c to show how thick the 
junction needs to be before sinusoidal behavior is restored (Ng = 20 is thick enough for sinusoidal 
behavior). The deviation is due to self-consistency. 

shape of this suppression is independent of the thickness of the junction, once the 
thickness is larger than a few bulk S coherence lengths. The exponential healing of 
F to its bulk value within the S occurs over the length scale of £s- Fitting the decay 
to an exponential gives a coherence length of £5 = 3.7a. Such a suppression of the 
superconductivity can be mimicked in the quasiclassical approach by introducing 
the so-called suppression parameter 75, but such calculations are rarely performed 
self-consistently. Here the decay and healing of the superconductivity is determined 
from the parameters of the microscopic Hamiltonian. It is important to note that 
for thin junctions L < 2£s, the value of F at the SN interface (plane 30) and the 
critical current of the junction depend strongly on the thickness of the junctions. 
Once L becomes large enough, the value of F near the SN interface is "frozen in" 
and the critical current is determined by the Thouless energy (for L on the order 
of £0) or by the normal metal coherence length £jy (for L larger than £/v). 
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Fig. 3. Phase deviation as a function of the plane number for a SNS junction with Nb = 20. The 
barrier begins at plane 31. Each curve corresponds to a different bulk phase gradient in the S leads, 
which increases from zero up to the critical current. Note how the largest phase gradient occurs 
at the center of the barrier. The total phase gradient (bulk plus deviation) is always positive. 

The next thing to notice is that the magnitude of F on the central plane 
of the barrier is determined solely from the (ballistic) N coherence length £jy = 
hvp I (2TrksT) « 35 — 40a (once the barrier is thick enough). As the barrier is made 
thicker and thicker, one can see how the pair-field decays exponentially with the 
thickness. For thinner junctions it is hard to detangle the two length scales, but for 
thicker junctions, one can clearly see the difference between £o and £jv in Figure ^. 

Finally, we can understand much about the current-phase relations, and the 
phase-deviation plot when we combine it with the proximity-effect plot of F. Note 
that as the magnitude of F decreases toward the center of the barrier, we find that 
the phase gradient is maximal at this central plane, since the magnitude of F is the 
smallest there. This is the only way that one can maintain current conservation from 
plane to plane. Hence, the critical current of the junction is found by the maximal 
phase gradient that can be sustained at the central plane of the barrier. This 
explains, in a simple way, why the (T > 0) critical current decays exponentially 




12 Optimizing the Speed of a Josephson Junction 



£ 0.06 




10 20 30 40 50 60 70 80 90 100 110 

Plane number 



Fig. 4. Anomalous pair field F versus plane number for a variety of junction thicknesses (1, 2, 4, 
6, 8, 10, 14, 20, 30, 40, and 60) and vanishing supercurrent. Note how the shape of F at the SN 
boundary does not depend on thickness once the junction is thicker than twice the bulk coherence 
length Ng > 8, and how the magnitude of F at the central plane is determined by the thickness 
of the junction and the N coherence length £jy 



(over the length scale £jv) as the thickness increases, since the magnitude of F 
decreases exponentially (for L large enough), and the maximal phase gradient is 
always the same order of magnitude. 

One final observation is in order. In the majority of the cases that we perform 
calculations for, there is little change in the proximity effect [i.e. in |-F(z)|] as one 
increases from zero supercurrent flow up to I c . In thick junctions, we see very 
little variation. As the junctions are made thinner, there are changes which are 
typically on the order of a few percent (but sometimes becoming large for very thin 
junctions). 

In addition to examining imaginary axis (equilibrium) thermodynamic proper- 
ties, we can also examine (equilibrium) dynamic properties. The most interesting 
quantity to calculate is the interacting density of states (DOS) at each plane in the 
self-consistent region. For a ballistic N junction, we know that the density of states 
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should show the presence of Andreev bound state peaks (which will be broadened by 
the coupling to the infinite leads and the sum over the transverse two-dimensional 
momentum), and we expect to see the density of states linearly approach zero at 
zero frequency within the barrier due to the presence of quasiparticles with vanish- 
ingly small longitudinal momentuma In Figure ||(a) we plot the DOS at the central 
plane of a moderately thick L = 10a SNS junction for zero supercurrent, moder- 
ate supercurrent (J = / c /4), and larger supercurrent (I = I c /2). As the current 
increases, we see that the peaks corresponding to the Andreev bound states move 
apart, due to the expected Doppler shift, since the Andreev bound states come in 
time-reversed pairs — one carries current to the right and one to the left. As super- 
current is passed through the junction, time-reversal symmetry is broken and the 
degeneracy of these states is lifted, because one state moves in the direction of the 
supercurrent flow, and one moves in the opposite direction. This can be clearly seen 
in the splitting of the peaks as supercurrent flows, which increases with increasing 
current. As the current is increased further, we find our computational algorithm 
breaks down on the real axis due to phase slips. In panel (b), we show the eveolution 
of the minigap (in an L = 5a SNS junction) as the attractive Coulomb interaction 
is reduced in magnitude. As expected for A/// scattering effects, we see that the 
minigap decreases in magnitude as the coupling strength decreases. In panel (c), 
we plot the local DOS as a function of plane number (with current / c /8 flowing 
through the junction) ranging from the center of the barrier (plane 34), to the SN 
interface (plane 30), to a length approximately £g inside the S (plane 25), to deep 
within the S (plane 10). Note how within the barrier, there is little dependence of 
the DOS on position (except for a small reduction near the minigap), and how the 
ABS leak far into the S (much farther than £g) and have relatively large spectral 
weight close to the SN boundary. Such effects can only be seen in self-consistent 
calculations. Finally in panel (d), we plot the current carrying DOS at the same 
small supercurrent as in (c) (/ = I c /8). The supercurrent is found by an integral 
of the current carrying DOS weighted by the appropriate Fermi factors. Note how 
the peaks of the current DOS correspond closely to the shifting Andreev peaks, 
confirming that the split peaks carry current in opposite directions. 

Surprisingly, the DOS does not decrease linearly to zer [see panels (a) and 
(b)], but rather shows the appearance of a minigap at the lowest energies. This 
behavior is different from what is expected from the quasiclassical situations, and 
we believe it arises from the fact that the superconductor has so short a coher- 
ence length that one cannot neglect A//i effects Jl3 (which are beyond quasiclassic 
methods) as illustrated most clearly in panel (b). As the junction is made thicker, 
the minigap decreases, and the linearly vanishing behavior is restored by the time 
L = 20a. The results for the DOS then agree with all of our expectations from the 
quasiclassical theory (once L is thick enough). This shows that our many-body for- 
malism is able to properly capture the regime described by the quasiclassical theory. 
But, as we will see below, it can go beyond this to account for correlation effects 
as well. Another point to emphasize is that the minigap depends only weakly on 
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the amount of current flowing through the junction. This is also contrary with the 
minigap that arises in quasiclassical situations, which show a decreasing minigap 
with supercurrent. This "theoretical spectroscopy" of a Josephson junction would 
be interesting to observe experimentally, and preliminary studies with tens of nm 
wide probes have already been performed.E3 It would be interesting to examine 
some of these junctions with the finer spatial resolution of an STM tip. 

The other dynamical property that yields much information about the junction 
is the normal state resistance R n . We calculate the normal state resistance by 
ignoring the superconductivity (setting F = 0), and using the Kubo formula for 
the conductivity. In performing our calculations, we neglect the effect of vertex 
corrections, which is exact in the large dimensional limit .Ej 

In a ballistic SNS junction, like the ones that we study here, the normal state 
resistance is essentially independent of the thickness L of the junction, and is given 
by the SharvintS resistance. Since the critical current depends exponentially on L 
for thick junctions (with characteristic length £;v)> we expect the I c R n product to 
decrease onca the thickness becomes larger than the N coherence length. Recent 
experimental showed anomalous behavior for SSmS junctions created out of Nb 
and heavily doped InAs. In particular, the characteristic voltage was seen to be 
much smaller than the Kulik-Omclyanchuk limit for a clean SNS junction. This 
should not come as a surprise, because there is a large Fermi velocity mismatch in 
these junctions, so the interfaces must have significant scattering even in "highly 
transparent" junctions. Whenever such scattering is present, one expects a reduc- 
tion from the Kulik-Omclyanchik limit. Here we examine the best possible scenario 
for the Kulik-Omclyanchuk limit, where there is no extrinsic scattering at the in- 
terface, and the only reductions to the I c R n product can arise from the proximity 
and "inverse" proximity effects or finite A/// scattering. 

In our junctions, we have no geometric constriction for the barrier, so the maxi- 
mal characteristic voltage, in the clean limit, is just the product of the bulk critical 
current with the Sharvin resistance of the infinite S leads. For our three-dimensional 
system, the product of these two is about one-half the Kulik-Omelyanchuk limit 
and equals 1.45A/e. What we find, is that even for junctions that have L < £jv ~ 
35 — 40a, the characteristic voltage is sharply reduced (for example, I c R n is an 
order of magnitude smaller than the thin-junction limit when L w 40a), similar 
to what is seen in experiment, but the effect isn't quite strong enough to produce 
the behavior seen in experiment (where the reduction in I c R n was estimated to be 
about two orders of magnitude) . The reason for the decrease in our calculations is a 
combination of the "inverse proximity effect" and the proximity effect, which reduce 
the superconductivity and thereby limit the amount of supercurrent that can flow 
through the junction. Such effects can only be found in self-consistent calculations. 
We will see in the next section how a charge redistribution at the SN interface (due 
to Fermi-level mismatch) reduces the characteristic voltage even further and is one 
possible way to explain the experimental behavior. 

Our work on SNS junctions has produced a number of interesting behaviors for 
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short-coherence-length Josephson junctions. We found that the critical current is 
determined by the maximal phase gradient that can be sustained over its "weakest 
link" — the central plane of the barrier. We also saw that because of the "inverse 
proximity effect" (where the superconductivity is reduced in the S as one approaches 
the SN interface) and the proximity effect (where superconductivity is reduced as 
one enters the N), the critical current has a strong dependence on the thickness 
of the junction. Since the normal-state resistance is essentially independent of the 
thickness, this provides a partial explanation for why the characteristic voltage of 
ballistic SNS junctions may be much smaller than predicted by analytic means. The 
reason is that the self-consistency, as seen through the proximity effects, sharply 
reduces the critical current, and thereby reduces the characteristic voltage. We 
verified, in a striking fashion, how the Andreev bound states split due to the presence 
of a supercurrent flowing through the junction, arising from the Doppler shift. 
Finally, we saw the appearance of a minigap for thin junctions, which gave way to 
the expected linear vanishing of the DOS as the junction is made thicker. 

4. SINIS junctions from double-barrier screened-dipole layers 

When two different metals are placed together in a hybrid structure, there is usually 
a mismatch of the Fermi levels. This means that one of the metals will have a 
chemical potential that differs from its equilibrium value. As a result charge will 
be redistributed at the interfaces between the metals until an overall equilibrium 
situation returns, where the charge deviates from the bulk values in the metals 
over a length scale (called the Debye screening length In) near the interfaces. This 
process should occur in SNS Josephson junctions as well. If the redistribution of 
charge is large enough, then it can create a SINIS junction, where the I layers are 
formed by the screened-dipole layers that appear at each SN interface. 

Hybrid SINIS junctions have been studied over the past two decades The 
basic idea is that one can combine the attractive properties of both SIS and SNS 
junctions by constructing hybrid SINIS junctions. The most heavily studied SINIS 
junctions are made out of the conventional Nb-Al-AlO^ process, with additional 
layers grown to create the more complicated structures. In some cases Nb is used 
in the barrier, which is a situation we will discuss in more detail in the next section. 
While much effort has been expended on trying to optimize junction properties, no 
one has been able to create reproducible SINIS junctions with high enough critical 
currents and characteristic voltages to be competitive with submicron, self-shunted 
tunnel junctions. These Nb-Al-AlO^-based junctions are described well by the 
quasiclassical approach,!!^ and we do not concern ourselves with them here. 

Instead, we concentrate on SINIS junctions where the I layers are self-consistently 
generated by screened-dipole layers in materials (such as the high-T c grain-boundary 
junctional or dopedE3 InAs junctions) where the Debye screening length is a few 
lattice spacings (rather than less than an Angstrom, which is seen in most conven- 
tional metals). In this regime, quasiclassical techniques are not applicable, because 
the redistributed charge profile requires one to take into account the modifications 
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of the Green's functions on length scales on the order of the Fermi lengthEj Xp- 
Quasiclassical techniques can only be attempted in materials where the screening 
length is so small that the potential from the screened-dipole layer can be rep- 
resented by a delta function. These spatially extended boundaries, however, are 
perfectly suited for our many-body approach. 

We use a combined quantum-classical technique to determine the self-consistent 
solution of the SINIS junctions. Quantum mechanics is used to determine the 
equilibrium charge density in each plane, subject to the given chemical potential 
and the current value of the local potential shift V.- L . Next, we determine the charge 
deviation 6n(z) = n{z) — neuik, where riBuik = tin or n$ if the plane is a normal 
metal or a superconductor, respectively. Once the charge deviations are known 
at each plane, we use classical electrostatics to determine the total potential shift 
at each plane, by summing the potential shifts (due to the constant electric fields 
emanating from each plane with a charge deviation) from every other plane. This 
then determines the total Coulomb potential at site i, Vf , and the potential energy 
shift Vi = eVf at each plane. Once the new potential energy shift is known, 
then we use the quantum-mechanical algorithm to calculate the charge densities at 
each plane, and repeat the process until the potentials and charge densities have 
converged to a fixed point. 

This procedure is time-consuming, and adds greatly to the computational effort 
needed to find a self-consistent solution. Fortunately, in most cases, the charge 
profile does not change much with either the supercurrent flowing through the 
system, or the temperature, so once the charge profile has been determined for one 
calculation, it serves as a good guess for the charge profile of the next calculation 
with similar parameters. Use of this "numerical annealing" strategy greatly reduces 
the number of computational cycles needed to achieve convergence of numerical 
results. 

In Figure 0, we plot the charge profile for Josephson junctions with L = 20a, 
a variety of Fermi- level mismatches AEp = Ep — Ep, and Id ~ 3a. The bulk 
superconductor filling is chosen always to be half filled ns = 1. In (a-b) the normal 
metal filling is also half-filled njv = 1, while in (c-d) it is much smaller tin = 0.01, to 
mimic the behavior in a doped semiconductor. Note how the charge deviation (and 
the corresponding potential barrier) grows as the Fermi-level mismatch increases 
in (a). We need only consider positive Fermi-level mismatch in this case, because 
particle-hole symmetry gives the same result for negative mismatches. One might 
note that the general shape of the charge profile looks the same for different mis- 
matches. This is confirmed in (b), where we rescale the charge profile, by dividing 
by AEp- The curves all collapse onto each other. We believe this occurs because 
the noninteracting density of states is quite flat near half filling for a simple cubic 
lattice. Indeed, when we look at the data in panel (d), which does not have any 
symmetry between positive and negative values, nor has the shape independent of 
the size of the Fermi level mismatch, we see that scaling does not hold for this 
data (which we believe is because the density of states is quite asymmetric and has 
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strong energy dependence near the Fermi level of the "doped semiconductor"). 

Recent Raman scattering experiments^ on Nb-InAs bilayers show a dramatic 
change in the Raman response as a function of the temperature as one cools through 
the superconducting transition of the Nb. One interpretation of these results, is 
that the effective Debye screening length, or equivalently, the charge profile at the 
interface, changes dramatically as one goes through T c . We see no evidence of this 
in our calculations, but we always fix the Debye screening length (or more precisely 
the optical dielectric constant eoo) in our calculations, and we see changes of the 
charge distribution no larger than a few percent as one cools through T c (the only 
exception is when AEp is on the order of T c or A, where sizable changes in the 
potentials, but not the screening lengths, are seen). It is our belief that if the spatial 
extent of the charge accumulation region is changing as one goes through T c , then it 
must require a proper treatment of the nonlocal screening effects, which is beyond 
our computational techniques. 

The next thing we examine is the proximity effect, shown in Figure |8[ In (a) we 
show the case with — 1 (AEp = 1,2,3,5, and 10) and in (b) we show tin = 0.01 
(AEp = 1, —1, —3, and —10). As the charge deviation increases, and the potential 
barrier increases, we see that the "inverse proximity effect" decreases as expected, 
and approaches the behavior expected for rigid boundary conditions (step function 
in A at the interface). But as the barrier increases further, due to larger mismatches, 
the "inverse proximity effect" gets larger and larger. This effect is due entirely to the 
modification of the local charge on the planes near the interface, and the fact that the 
Debye screening length is chosen to be a few lattice spacings. Note how the presence 
of a scattering charge barrier at the interface sharply reduces the superconductivity 
within the N. This is a consequence of the reduced transparency due to the charge- 
redistribution-induced scattering. In the case of the "doped semiconductor" , we 
find that for small mismatches |A.Ef| < 1, the proximity effect does reproduce the 
expected rigid boundary conditions, with the added feature of small oscillations 
(due to Fermi-length effects) near the interface. But as the mismatch is made more 
negative, one reproduces the same kind of enhanced "inverse proximity effect," as 
was seen at half filling, due to the charge redistribution within the superconductor. 
Note that the step-function boundary conditions for small mismatch are most likely 
due to the fact that there is a huge Fermi-wavevector mismatch in this junction, 
which causes large scattering, and reduces the transparency of the junction to the 
point where the conditions for rigid boundary conditions hold. The appearance of 
the small oscillations would not be seen in a quasiclassical approach, that averages 
over short length scales. 

Our final result is the characteristic voltage versus the junction thickness for 
three different Fermi level mismatches in Figure ^. The critical current decreases 
monotonically with the thickness, in nearly all cases. The normal state resistance, 
can be enhanced dramatically for intermediate sized junctions, where the thickness 
is less than twice the Debye screening length. Once the thickness is greater than 
21 d, the resistance is virtually independent of the thickness, as expected for ballistic 
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junctions. It is mainly this anomalous behavior at intermediate thicknesses, which 
leads to the nonmonotonic behavior of the characteristic voltage. When we compare 
I c R n to the SNS case, we see that for thicker junctions, the characteristic voltage is 
sharply reduced when there is additional scattering due to a screened-dipole layer. 
This arises mainly from the fact that the critical current is reduced more rapidly 
with thickness for thicker junctions. But, surprisingly, there is an intermediate 
regime, of relatively thin junctions (similar to grain-boundary junctions in high- 
T c ), where the characteristic voltage can be enhanced by the charge redistribution 
(AE F = 5 and N B = 3). 

We find that when a screened dipole layer develops at a SN interface due to 
Fermi-level mismatch, it usually produces a sharp reduction in the characteristic 
voltage, because the critical current is reduced much more than the normal-state 
resistance (because of the large reduction of the proximity effect due to the scatter- 
ing from the extended charge barrier). This is clearly not advantageous for digital 
electronics applications (or for power applications in multigrain high-T c tapes). It 
appears that junction qualities should be improved if the Fermi-level mismatch can 
be reduced by an educated choice of materials used in the junctions. 

5. SNSNS junctions with enhanced I c R n products 

We saw in the previous section that if we allow charge to redistribute due to a 
Fermi level mismatch, then it usually has a deleterious effect on the characteristic 
voltage. If we recall that the critical current (of a SNS junction) is determined by 
the maximal phase gradient that the weakest plane can sustain, and that the normal 
state resistance is dominated by the Sharvin resistance in clean systems, then one 
way to enhance the characteristic voltage is to improve the superconductivity within 
the central planes of the barrier, while maintaining quantum coherence throughout 
the whole junction. This motivates the consideration of SNSNS junctions where 
one replaces the central planes of the N barrier with S. Recent experimental have 
indicated that one can get a d ramatic rise in the characteristic voltage, which has 
been confirmed theoretically^!. 

Here we will explore a simple system, where we take the 20 barrier planes and 
replace the middle 6 planes with S. In Figure [l^, we plot the proximity effect within 
this junction, and see that the proximity effect is enhanced due to the extra super- 
conducting planes (the minimal F is about twice as large in the SNSNS junction 
than in the SNS junction). Because of this, the critical current is sharply enhanced, 
increasing by a factor of two over that of a similar SNS junction with L — 20a. Since 
there are additional interfaces, which will add to the scattering, the resistance will 
be increased too, but since it is dominated by the Sharvin contribution, it remains 
essentially the same. The net effect is a near doubling of the characteristic voltage 
to 0.775A/e (I c R n = 0.46lA/e for the L = 20a SNS junction); it still remains much 
lower than the maximal analytic limit of nA/e. 

Even more interesting is a plot of the many-body DOS as a function of position 
within the SNSNS plane. The leftmost SN boundary is at plane 30, the inner NS 
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boundary is at plane 37, the inner SN boundary at plane 44, and the rightmost NS 
boundary is at plane 50. One can see interesting structure in the Andreev bound 
states. The main peak at approximately ±0.6A remains at about the same energy in 
the interior S planes and in the N. As we move into the N, away from the innermost 
S planes, there is a small spectral weight shift to lower energies, near the minigap 
in the DOS, and at about 0.9A where a small peak develops in the N. Note that 
even though the L = 20a SNS junction has no minigap, we find a sizable minigap 
here, because it is determined by the A/// scattering at each of the (now four) SN 
interfaces. As we move deep into the superconductor (approximately 5£s from the 
bulk SN interface at plane 10), we still see a remnant of the Andreev bound state 
peak in the superconductor at about 0.6A, and a recovery of the expected BCS 
density of states elsewhere. Similar problems have also been investigated by Lodder 
and coworkers .lj 

We see that there are complicated structures that can lead to dramatic enhance- 
ments of the characteristic voltage. In particular, replacing some of the central 
planes of the N barrier by S will enhance the proximity effect and thereby enhance 
J c , and will also typically enhance the normal state resistance due to extra inter- 
faces. This result has been seen already in experiment.E3 What is interesting, is 
that if the S layers are thin enough, then the general character of the ABS remain 
unchanged within the hybrid junction, and the net effect is an increase in the char- 
acteristic voltage of the junction. What may prove to be more problematic in these 
junctions is whether or not they are stable to changes in temperature, which is 
necessary for implementation within digital electronics. 

6. SCmS junctions tuned through a metal-insulator transition 

The final type of junction we will consider here is quite different from the others, and 
has not been analyzed with quasiclassical approaches, because it lies well outside the 
region of validity for quasiclassics. It is a SCmS junction, where the correlated metal 
is described by the Falicov-Kimball model, tuned to lie close to the metal-insulator 
transition in the bulk. Depending on precisely where one sits in the phase diagram, 
one can tune the metal-insulator transition to occur either at a fixed thickness, by 
increasing the interaction str eng th, or at a fixed interaction strength, by increasing 
the thickness L of the barrierEil. 

The Falicov-Kimball model is always a non-Fermi liquid, because the presence 
of static disorder always creates a finite lifetime to putative quasiparticles at the 
Fermi surface. c3 Initially, this behaves like a "disordered" Fermi liquid, but as the 
interaction increases in strength, there is a transition in the many-body DOS to 
a pseudogap structure, where there is a dip at the Fermi level, followed shortly 
thereafter by a true gap developing at the chemical potential for Ufk > 4.9. If 
Ufk > 4.9, then a metal-insulator transition occurs as a function of the thickness 
of the barrier (since a single plane barrier [Nb = 1] is never insulating for transport 
perpendicular to the plane). Similarly, if we fix the barrier thickness, then there is 
a metal-insulator transition as a function of Ufk- 
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We first examine the proximity effect for the case where Nb = 5, which is just 
slightly larger than the bulk coherence length of the S. When Ufk is small, the 
junction behaves similarly to a SNS junction, with both "inverse" and conventional 
proximity effects. As Ufk increases to larger than 6, we first see the "inverse 
proximity effect" reduced, with the development of small oscillations due to "Fermi- 
surface effects," but the most striking development is of sharp oscillations, with a 
large amplitude within the barrier that are tied to the SB interfaces. We believe 
that these oscillations arise from properties of an insulating barrier being brought 
in contact with a superconductor. Oscillations develop at each SB interface over a 
length scale on the order of £oj which can interfere if the barrier is thin enough. We 
find that as the barrier is made thicker, the oscillations quickly die off as a function 
of distance from the SB interface (but the amplitude near the interface remains 
large). 

As disorder scattering isadded to the barrier of a junction, a new energy scale, 
called the Thouless energjH £>rh starts to play a role in determining properties 
of the junction. The Thouless energy is defined to be E^t = TiV/L 2 , with T> the 
classical diffusion constant. One of the most important predictions of the quasi- 
classical approach, is that the interacting DOS develops a minigap on the order 
of i/rh near the chemical potential for a "disordered" Fermi liquid in a confined 
mesoscopic geometry!!! This gap increases in size as the disorder increases until it 
reaches a maximum at some fraction of the superconducting gap, and then the dirty 
limit takes over. This modification of the DOS in the normal metal, from essen- 
tially constant, to possessing a low-energy "minigap" is one of the most interesting 
consequences of the proximity effect, when analyzed in a quasiclassical approach 
(especially its evolution with increasing disorderEj) . 

The behavior in a Cm is quite different. We plot in Figure [l3] the many-body 
DOS at the central plane of a L = 10a SCmS junction (i.e., for the plane at the 
center of the barrier), for values of Ufk ranging from 0.1 to 2. When Ufk = 0, 
a minigap appears in the DOS as seen in Figure |^. As scattering is introduced, 
the minigap shrinks, eventually disappearing and becoming a pseudogap as Ufk 
increases. For small Ufk we see the characteristic development of a minigap at 
low energy, but we see another energy scale enters as well at about five times the 
minigap, where a "soft" pseudogap appears. As U fk increases further, this picture 
evolves, with the gaps remaining intact, but narrowing, until a critical value of 
Ufk ~ 0.8 is reached, and the minigap disappears. The DOS continues to fill in 
at low energies, with the remnants of a pseudogap remaining until Ufk becomes 
large enough that the DOS loses all low-energy structure. Note further, that as the 
scattering increases, by increasing Ufk, we find that the Andreev bound states, 
occurring just below the superconducting gap, first broaden dramatically, then get 
washed out of any structure as the correlations increase further. In fact, it is likely 
that the destruction of the minigap is occurring because of the extensive broadening 
of the Andreev bound states. 

What is interesting, is that all of the behavior of the closing of the minigap occurs 
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well before the pseudogap develops in the bulk many-body DOS of the Falicov- 
Kimball model. Another surprise is that if we extract the diffusion constant from 
the normal-state resistance we find that the size of the gap is much smaller than 
the quasiclassical predictions of 3.12i?Th- In addition, we found that the size of the 
gap does not scale inversely as L 2 , as it must in the quasiclassical theory. Hence the 
proximity effect in a correlated metal has very different characteristics than what 
is seen in a disordered Fermi liquid. 

We now turn to the question of optimizing the characteristic voltage of a SCmS 
junction. In Figure [li], we show plots of the characteristic voltage versus Ufk for 
four different thicknesses of the barrier L = a, 2a, 5a, and 20a. The dashed line is 
the Ambegaokar-Baratoff predictionH. A thin junction (L = a) behaves pretty much 
as expected. I c R n lies below the planar junction limit for small Ufk and increases 
modestly as the scattering increases, reaching an optimization near Ufk — 2. As 
one goes through the metal-insulator transition, the characteristic voltage saturates, 
and becomes independent of Ufk ■ It lies about 15% below the Ambegaokar-Baratoff 
prediction though, because of the "inverse proximity effect" and the averaging over 
an anisotropic "Fermi surface." When we double the number of planes to two, we 
find similar behavior at small Ufk with a mild optimization of I c R n near Ufk = 3, 
but the characteristic voltage continues to increase without any apparent bound as 
Ufk increases in the insulating regime. This is quite different from the Ambegaokar- 
Baratoff prediction. Of course, we expect the junctions to be hysteretic when the 
system is a good insulator. The behavior for the intermediate-sized junction L = 5a 
is the most interesting. In the metallic region, the characteristic voltage decreases 
monotonically with Ufk until we hit the metal-insulator transition. At that point 
I c R n increases by more than a factor of 2, shows a mild optimization near Ufk = 7, 
and then decreases for larger Ufk- This illustrates the proposition that I c R n can be 
optimized near the metal-insulator transition, and be higher than the Ambegaokar- 
Baratoff prediction. We examine the final case of a thick junction L = 20a in panel 
(d). Here, the characteristic voltage decreases sharply with U fk, with the rate being 
fastest near the metal-insulator transition. We believe the very low values on the 
insulating side occur because the normal state resistance, being thermally activated, 
has strong temperature dependence even at these low temperatures, which leads to 
the extremely low values of the characteristic voltage. 

Finally, we produce a "Thouless plot" of our data, by fixing Ufk and study- 
ing the characteristic voltage as a function of the barrier thickness. We consider 
three cases: (i) Ufk — 2, a strongly disordered metal; (ii) Ufk = 4, a pseudogap 
metal; and (nil Ufk = 6, a correlated insulator with a small gap. A quasiclassi- 
cal predictionE3, shows that the characteristic voltage smoothly changes from the 
Kulik-Omelyanchuk limfola to being proportional to -&rh for thick junctions. The 
crossover is a universal curve in the ratio of the Thouless energy to the superconduct- 
ing gap. We plot this behavior in Figure [IS The Thouless energy is determined by 
first fitting the exponential decay of the supercurrent with L to extract the barrier 
coherence length fg. The Thouless energy is then E^h = ^irksT^/ L 2 . 
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Note that both metallic curves (Ufk — 2 and 4) fall on essentially the same 
curve, illustrating the "universality" of the quasiclassical prediction.E3 This may 
be a little surprising, because the pseudogap metal is clearly outside of the realm 
of the quasiclassical approach. The correlated insulator has a nonuniversal shape, 
that deviates sharply from that of the metals. The characteristic voltage is initially 
quite flat for thin junctions, and then is reduced sharply with thickness for thicker 
junctions. Because of this behavior, we expect the Josephson coupling to be very 
sensitive to thickness when L is larger than the crossover between the metal and 
the insulator (where the curve decreases sharply). In this region, the Josephson 
coupling can be much larger for slightly thinner barrier regions, and much smaller 
for slightly thicker regions. This behavior can look like pinholes, but will occur, even 
if the barrier is homogeneous, and has a small variation in its thickness (perhaps 
even as small as one monolayer). We believe this kind of behavior may occur in 
high-T c junctions, which have much larger spreads in their junction parameters even 
in well controlled fabrications runs. 

We found a number of new behaviors in correlated-barrier Josephson junctions. 
In particular, there is an anomalous enhancement of the proximity effect at the SB 
interface, as the barrier is made more insulating. In addition, we found a wide vari- 
ety of different characteristic voltages, with an optimization on the insulating side 
of the metal-insulator transition for intermediate-thickness junctions. We also saw 
remarkable behavior in the interacting DOS, where the quasiclassical prediction of 
a minigap, disappears as correlations increase, and broaden the ABS in the DOS; 
first the minigap is replaced by a pseudogap, and then it becomes a structureless 
flat background as the correlations increase further; this effect is most likely due 
to "lifetime" effects that broaden the Andreev bound states into wide, featureless 
resonances. Finally, we found that correlated insulating barriers can have an in- 
trinsic pinhole effect, when the thickness lies close to the critical thickness of the 
metal-insulator transition. 

7. Conclusions 

We have illustrated many of the novel results that can be analyzed with a new ap- 
proach toward modeling Josephson junctions, that are relevant for short-coherence- 
length junctions with either correlated barriers, or barriers with potentials that have 
features on atomic length scales. In this review, we have emph a sized the physical 
properties, since the formalism has been developed elsewhere. We examined 
many different types of junctions, ranging from SNS, to SNSNS, to SINIS, to SIS- 
mlS, to SCmS. We found a number of interesting results that are not seen in the 
quasiclassical theory. These include (i) a sharp reduction in the characteristic volt- 
age when a mismatch of Fermi levels causes a screened-dipole layer; (ii) a large 
enhancement of I c R n for SNSNS junctions; (iii) the disappearance of the minigap 
due to increased scattering in correlated barriers; (iv) an optimization of I c R n for a 
correlated insulator close to the metal-insulator transition; and (v) the prediction of 
an "intrinsic pinhole" effect for correlated insulating barriers due to deviations from 
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universality in the Thouless plot. These results can explain a number of anoma- 
lous results seen in experiments on Josephson junctions ranging from ballistic SSmS 
junctions to high-T c -based junctions. 

The use of dynamical mean field theory in modeling properties of inhomogeneous 
structures like Josephson junctions allows us to calculate properties of a number 
of systems that lie outside of the successful quasiclassical approach which include 
short-range potentials, correlated systems, and structures with thicknesses on the 
order of the Fermi length. These techniques can easily be used for other types of 
devices, such as hybrid superconducting and ferromagnetic structures, with poten- 
tial applications to spintronics and quantum computing. In future work, we will 
generalize these results to allow nonequilibrium properties (such as the I — V char- 
acteristic) to be calculated, to allow unconventional pairing symmetry (like d-wave) 
for the superconductors for a better simulation of high-T c devices, and to allow 
ferromagnetic states for simulating spintronic devices. 
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Fig. 5. Many-body density of states as a function of frequency for L = 10a. (a) The solid 
line is with no supercurrent, the dotted line with moderate supercurrent Vc/4, and the dashed 
line for larger supercurrent / c /2, at the central barrier plane. Note how the DOS has a minigap 
at small frequency, and how the peaks, corresponding to Andrecv bound states, move apart as 
supercurrent flows through the junction, (b) Detail of the minigap for U = —2 (solid line) and 
U = —1.7 (chain-dotted line) and L = 5a. Note how the reduction of the size of the minigap for 
weaker interactions indicates that the likely source of the minigap is the so-called A/fi scattering, 
(c) Local DOS as a function of plane number for I = I c /S. Note how the ABS leak far into 
the superconductor (distances much larger than £g), a result that could not be seen with non- 
sclf-consistent calculations, (d) Current-carrying density of states at the central barrier plane for 
/ = / c /8. 



Optimizing the Speed of a Josephson Junction 27 



-I I I I I I I I I | I I I I I I I I I | I I I I I I I I I | I I I I I I I I I | I I I I I I I I I | I I I I I I I I I | I- 




10 20 30 40 50 60 

Junction Thickness [a] 

Fig. 6. Scmilogarithmic plot of the characteristic voltage versus junction thickness L for ballistic 
SNS junctions. The dashed line is the clean planar junction limit. Note how the characteristic 
voltage is suppressed dramatically below the clean planar limit due to the "inverse proximity 
effect" and the proximity effect. 
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Fig. 7. Charge redistribution in a SINIS Josephson junction with (a-b) tin = 1 and (c-d) 
"at = 0.01. The charge deviation 5n(z) is plotted for a number of different values of the Fermi 
level mismatch AEp (AEp = 0.1, 0.5, 1, 2, 3, and 5 for n N = 1 and AEp = 1, -1, -3, -10 for 
"at = 0.01). The values of AEp decrease from the bottom line to top (at plane 35) in panels (a, 
c, and d). Panels (b) and (d) show rescaled plots, where we divide the charge deviation by AEp. 
Note how the charge redistribution scales for the half-filled (metallic) case, but does not for the 
semiconductor case. The distributions are symmetric for positive and negative AEp in (a-b), but 
not in (c-d). 
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Fig. 8. Proximity effect for (a) njv = 1 and (b) njv = 0.01. Note how, in (a), the "inverse 
proximity effect" first is reduced as the barrier height increases, but then the "inverse proximity 
effect" gets larger, as the spatially extended barrier increases in magnitude. In (b), we recover rigid 
boundary conditions for small mismatch (plus Fermi- length oscillations), which then disappear as 
the barrier gets larger. Note also how the proximity effect is sharply reduced as the barrier gets 
larger. 
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Fig. 9. Scmilogarithmic plot of the characteristic voltage versus junction thickness L for ballistic 
SINIS junctions. The dashed line is the clean planar junction limit. Note how the characteristic 
voltage is suppressed dramatically below the clean planar limit for thicker junctions, but can be 
enhanced when the junction thickness is on the order of the Debye screening length (AEp = 5 
and Ng = 3). Note that a small barrier AEp = 1 does not change much from the case with no 
Fermi-level mismatch. 
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Fig. 10. Proximity effect for a SNSNS junction (solid line), where the central 6 planes of a SNS 
junction (with Nb = 20) are replaced by S planes compared to the proximity effect in the SNS 
junction (dashed line). Note how the proximity effect is sharply enhanced within the central S part 
of the junction due to the enhanced superconductivity in that region. This enhancement leads to 
a significantly enhanced I c (in fact, the ratio of the enhanced I c Rn value of about 1.7 is close to 
the ratio of the minimal value of F in the two junctions, which is about 1.8). 
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Fig. 11. Many body density of states as a function of position in the SNSNS junction. The solid 
curve is at plane 40, at the center of the innermost S layers; the dashed curve is at plane 35 just 
inside the N metal; the chain-dotted curve is at plane 31, at the SN boundary, and the dotted line 
is at plane 10, approximately 5£s away from the SN boundary within the bulk S. Note how the 
main Andreev peak at 0.6A remains essentially unchanged within the interior of the junction (and 
even has some weight deep within the S). We see an additional development of spectral weight 
near both the minigap and near the bulk gap in the N relative to the innermost S layers. The 
appearance of the minigap for this junction arises from the increased A/fj, scattering due to the 
additional two SN (NS) interfaces. 
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Fig. 12. Proximity effect for a SCmS junction with L = 5a. Note how the shape of the anomalous 
average smoothly changes from that of a SNS junction, with a large "inverse proximity effect," to 
that of a SI interface, with a relatively fiat F in the superconductor (plus small oscillations). The 
interesting result is the large amplitude oscillation occurring within the barrier, which dies off in 
a length scale on the order of £o- The curves alternate from solid and dashed lines for Upn = 1, 
2, 3, 4, 5, 6, 7, 8, 10. 
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Fig. 13. Many-body density of states at the central plane of a L = 10a SCmS junction. The 
correlation strength Uppc is increased from 0.1 to 2. Note how we first see a minigap that disappears 
as the correlations increase and becomes a pseudogap, before vanishing entirely into a structureless, 
fiat "background." 
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Fig. 14. Semilogarithmic plot of the characteristic voltage versus Ufk for four different junc- 
tion thicknesses: (a) L = a; (b) L = 2a; (c) L = 5a; and (d) L = 20a. The dashed line is 
the Ambegaokar-Baratoff prediction of I c R n = 1.57A/e. Note how in (a) we do reproduce the 
Ambegaokar-Baratoff prediction of the characteristic voltage being independent of Ufk in the 
insulating regime, but the value is somewhat smaller in magnitude. In (b), we see that I c R n 
appears to increase without limit. In (c), we find an optimization, above the Ambegaokar-Baratoff 
prediction, just on the insulating side of the metal-insulator transition. And in (d), we find I c R n 
decreases dramatically with Ufk- 
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Fig. 15. Thouless plot for Ufk = 2, 4, and 6. The characteristic voltage is plotted against the 
Thouless energy. Note that both metallic cases (Upx = 2 and— 4) do fall on a nearly "universal" 
curve, which is similar in shape to that seen in quasiclassicsEHI The insulating curve, however, 
deviates from this form, having a flat curve for thin junctions, and then a sharp decrease after the 
thickness-tuned metal-insulator transition. This behavior leads to the possibility of what we call 
an intrinsic pinhole effect. 



